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ABSTRACT 

We study the spectra of all long Gamma Ray Bursts (GRBs) of known redshift detected by the Fermi satellite. Their fluxes and 
fluences are large enough to allow a time dependent study of their spectral characteristics in the 8 keV-I MeV energy range. We 
find that the peak energy iipeak of their EL(E) spectrum correlates with the luminosity in a remarkable tight way within individual 
bursts. This time resolved iipeak - Uso correlation is very similar for all the considered bursts, and has a slope and normalisation 
similar to the analogous iipeak - hso correlation defined by the time integrated spectra of different bursts detected by several different 
satellites. For a few of the considered GRBs, we could also study the behaviour of the iipcak - i^iso correlation during the rising and 
decaying phases of individual pulses within each burst, finding no differences. Our results indicate the presence of a similar physical 
mechanism, operating for the duration of different GRBs, linking tightly the burst luminosity with the peak energy of the spectrum 
emitted at different times. Such a physical mechanism is the same during the rise and decay phase of individual pulses composing a 
GRB. These results, while calling for a robust physical interpretation, strongly indicate that the iipeak - ^iso spectral energy correlation 
found considering the time integrated spectra of different bursts is real, and not the result of instrumental selection effects. 
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' 1. Introduction 

One of the key properties of the prompt emission of Gamma 
Ray Bursts (GRBs) that is still poorly understood concerns the 
spectral-energy correlations found considering the time inte- 
grated spectra of bursts of known redshift, and for which we 
can define the peak energy of the spectrum, fipeak, in a EL(E) 
, representation. /Speak correlates with the isotropic luminosity Liso 
' (Yonetoku et al. 2004), with the isotropic energy Eiso (Amati et 
" al. 2002) and more tightly with the coUimation-corrected en- 
ergy Ey (Ghirlanda, Ghisellini & Lazzati 2004). There are two 
very strong motivations to study these correlations: one is to un- 
derstand their physical origin, that can disclose a still not un- 
derstood basic property of GRBs (Yamazaki, loka & Nakamura 
2004; Lamb, Donaghy & Graziani 2005; Rees & Meszaros 2005; 
Levinson & Eichler 2005; Toma et al. 2005; Eichler & Levinson 
2006, 2004; BarbielHni et al. 2006; Thompson 2006; Ryde et al. 
2006; Giannios & Spruit 2007; Thompson, Meszaros & Rees, 
2007; Guida et al. 2008; Panaitescu 2009) and the other is the 
possibility to use these correlations to standardise the GRB ener- 
getics, making them cosmological tools (Ghirlanda et al. 2004a, 
2006, 2006a; Firmani et al. 2005, 2006, 2007; Xu, Dai & Liang 
2005; Liang & Zhang 2005, 2006; Wang & Dai 2006; Qi, Wang 
& Lu 2008; Li et al. 2008; Liang et al. 2008). 

The debate about the reality of these correlations is hot, 
since some authors pointed out that they can be the result of 
observational selection effects (Nakar & Piran 2005; Band & 
Preece 2005; Butler et al. 2007, Buder, Kocevski & Bloom 2009; 
Shahmoradi & Nemiroff 2009) while others argue that selection 
effects, even if surely present, play a marginal role (Ghirlanda et 
al. 2005, Bosnjak et al. 2008, Ghirlanda et al. 2008; Nava et al., 
2008; Krimm et al. 2009; Amati et al. 2009). 



One possibility to get some insight on this issue is to study 
individual, bright bursts to see if, during the prompt phase, the 
luminosity and peak energy at different times correlate. If they 
do, and furthermore if the slope of this time resolved correlation 
(indicated E' , - L' hereafter) is similar to the time integrated 

' peak ISO ^ ^ 

£^peak - Liso coiTclation found among different bursts, then we 
should conclude that the spectral energy correlations are surely 
a manifestation of the physics of GRBs, and are not the result of 
instrumental selection effects. 

Some attempts have already been done. Liang, Dai & Wu 
(2004) considered BATSE bursts without known redshifts and 
showed the presence of a correlation between the (observer 
frame) peak energy and the flux within individual bursts which 
they interpret as suggestive of a physical origin of the /Speak - Liso 
coiTelation holding among the GRBs with measured redshift. 
However, to compare direcly the E'^^^y^ - L[^^ correlations of indi- 
vidual GRBs with the fipeak - i'iso correlation defined with time 
integrated spectra it is necessary to know the redshift (which is 
instead unknown for most of BATSE bursts). Recently, Firmani 
et al. (2009) considered Swift bursts of known redshift finding 
a rather strong E' , - V. correlation within individual GRBs. 

^ peak ISO 

Having the redshift, they could directly compare the time re- 
solved correlation of different bursts, finding that the ensemble 
of data points in the iipeak - i^iso plane shows a correlation similar 
to that defined with the time integrated spectra of the same burst 
sample. The Burst Alert Telescope (BAT) onboard Swift, how- 
ever, with its limited energy range (15-150 keV), is not particu- 
larly suited for GRB spectral analysis, especially when dealing 
with time resolved spectra. To overcome this limitation, Ohno et 
al. (2009) combined the Swift-QM and Suzaku-WAM spectral 
data to study the spectral evolution of GRB 061007 and inves- 
tigate the time evolution of the iipeak - i'iso correlation within 
the two pulses of this burst. They found that the time resolved 
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Table 1. Fermi long GRBs with redshift. We mark in boldface the GRBs with a detection also in the LAT instrument on-board 
Fermi* The energy peak flux of GRB 080916C in the 20keV-10MeV range is from Golenetskii et al. 2008 



GRB 


z 


a 


^peak 




P 


F-6 


range 


GCN 


^180.52 


£150,52 








keV 




ph/s/cm^ 


erg/cm- 


keV 


number 


erg/s 


erg 


080810(549) 


3.35 


-0.91±0.12 


313.5+3.6 




1.85+0.16 


6.9+0.5 


50-300 


8100 


7.84 


33.2 


080905(705) 


2.374 


-1.75±0.12 






0.21+0.02 


0.04+0.003 


20-1000 


8205 






080916(406) 


0.689 


-0.9±0.1 


109+9 




4.5+0.7 


15+5 


25-1000 


8263 


0.142 


2.25 


080916(009) 


4.35 


-0.91+0.02 


424+24 


-2.08+0.06 


1.2e-5* 


190 


8-30000 


8278 


190 


563 


080928(628) 


1.692 


-1.80+0.08 








1.5+0.1 


50-300 


8316 






081007 


0.529 


-1.4±0.4 


40+10 




2.2+0.2 


1.2+0.1 


25-900 


8369 


0.041 


0.172 


081222(204) 


2.77 


-0.55+0.07 


134+9 


-2.10+0.06 


14.8+1.4 


13.5+0.8 


8-1000 


8715 


20.6 


35.4 


090323(002) 


3.57 


-0.89+0.03 


697+51 




12.3+0.4 


100+1 


8-1000 


9035 


47.2 


338 


090328(401) 


0.736 


-0.93+0.02 


653+45 


-2.2+0.1 


18.5+0.5 


80.9+1 


8-1000 


9057 


1.96 


21.2 


090423(330) 


8.2 


-0.77+0.35 


82+15 




3.3+0.5 


1.1+0.3 


8-1000 


9229 


18.8 


10.2 


090424(592) 


0.544 


-0.90+0.02 


177+3 


-2.9+0.1 


137+5 


52+1 


8-1000 


9230 


2.12 


4.48 


090618(353) 


0.54 


-1.26+0.04 


155.5+11 


-2.5+0.25 


73.4+2.0 


270+6 


8-1000 


9535 


1.0 


25.7 



pulses also satisfy the ZSpeak - ^-iso correlation defined by time 
integrated spectra. A more systematic analysis of the time re- 
solved spectral properties of Swift-Suzaku GRBs (Krimm et al. 
2009) shows that individual pulses within a GRB are consistent 
with the fipeak - £^iso Correlation defined by the time averaged 
spectra. They consider the spectra integrated over the duration 
of individual pulses. In this case /Siso is computed on different 
integration timescales. They find that the pulses follow a corre- 
lation parallel to the fipeak - £^iso defined with time integrated 
spectra, but with a higher normalization. Instead, the compar- 
ison of time resolved spectra with the /Speak - ^iso correlation 
is independent of the spectral integration time since this corre- 
lation involves the luminosity Liso rather than the energy Eiso- 
It is therefore important to study the presence of a fipeak - i'iso 
correlation (as done by Ohno et al. 2009 for a single event) by 
concentrating on GRBs with measured redshift, in oder to com- 
pare their time resolved E' , - L' correlation with that defined 

^ peak ISO 

with the time integrated spectra. In particular we aim at studying 
how single GRBs evolve in the fipeak -i'iso plane rather than con- 
sidering them globally (as in Firmani et al. 2009). Furthermore, 
we would like to study the rise and decay phases of individual 
pulses. 

The Gamma ray Burst Monitor (GBM, Meegan et al. 2009) 
onboard the Fermi satellite covers a wide energy range (8 keV - 
30 MeV) and, although slightly less sensitive than 5w//f/BAT, 
it is better for studying the spectral properties of the prompt 
emission of GRBs. In addition, the Large Area Telescope (LAT) 
sensitive in the 0.1-100 GeV energy range can complement the 
spectral information for the few bursts it can detect. In one year 
of operation (up to the end of July 2009), Fermi/GBM detected 
about 200 bursts and for about half of them the (time integrated) 
spectral analysis returned a well defined E°^\ . 

^ J peak 

In this paper we study Fermi GRBs, selecting those of known 
redshift, to be able to compare their different evolutionary tracks 
in the fipeak - i'iso plane with the correlation defined by the time 
integrated spectra (see e.g. Ghirlanda et al. 2009 for a recent 
compilation of the iSpeak - ^-iso correlation). For the brightest 
bursts we study if the rising and decaying phases of individual 
pulses behave differently in the /Speak - i-iso plane, since this can 
give important clues for our physical understanding of the emis- 
sion mechanism operating during the GRB prompt phase. 

The paper is organised as follows: in §2 we describe our 
Fermi GRB sample, whose time-integrated spectral properties 
are presented in §3 and compared to the fipeak - £^iso and fipeak - 
Liso correlation defined by pie-Fermi GRBs. In §4 we describe 
the time resolved spectral analysis whose results are given in §5. 



In §6 we discuss our findings and we draw our conclusions. A 
standard cosmology for a flat universe with /iq = - 0.7 is 
assumed. 

2. The sample 

We considered the GRBs detected by the GBM (up to the end 
of July 2009) of known redshift. They are 13 events. Among 
these GRB 090510 (z=0.903, Rau et al. 2009) is a short burst, 
having an observer frame duration of less than 2 s, and will not 
be considered here. 

Tab.[T]lists the 12 long GRBs, their time integrated spectral 
parameters (Col. 3 to Col. 8) and the derived isotropic luminos- 
ity (Liso, Col. 10) and isotropic energy (Liso, Col. 1 1), computed 
in the rest frame 1 keV - 10 MeV energy range. The spectral 
parameters in Tab.[T]have been collected from the literature (ref- 
erences are given Col. 9): they were obtained through the anal- 
ysis of the time integrated spectrum extracted from the GBM 
data. In two cases (GRB 080905 - Bhat et al. 2008; and GRB 
080928 - Paciesas et al., 2008) the time integrated spectrum is 
fitted by a single power law and, therefore, the peak energy is 
unconstrained. In five cases the time integrated spectrum is mod- 
eled with a power law ending with an exponential cutoff at high 
energies. In the remaining five cases it is modeled by a Band 
function. 

Two of the three GRBs detected by the LAT (in boldface 
in Tab. [TJ belong to the latter group (i.e GRB 080916C, Tajima 
et al., 2008; Abdo et al. 2009; and GRB 090328, Cutini et al., 
2009). These bursts show a high energy power law component 
and their observed peak energies Lp^^k '■^^ largest of the 
sample. Note that Tab. [T] lists also the most distant burst: GRB 
090423, with z=8.2 (Tanvir et al., 2009). 

3. Time integrated spectra: the Epsak - £^iso and the 

£^peak - ^iso correlations 

First, we check the consistency of Fermi bursts with the Lpeak - 
Liso and the Lpeak - Liso correlations defined by the time inte- 
grated spectra of GRBs detected by other instruments. The most 
updated pis-Fermi sample of GRBs with known z and L°^^j, con- 
tains 100 objects detected by diff'erent instruments. Fig.^shows 
these bursts (grey filled circles) in the Lpeak - Liso and Lpeak - Liso 
planes. In both planes they define a strong correlation (the prob- 
ability that the correlation is by chance is reported in Col. 4 of 
Table|2]l, confirming recent analysis (Nava et al. 2008; Ghirlanda 
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et al. 2009). The slope and normalization of the fit of these corre- 
lations with a power law are reported in Table |2] To show where 
the Fermi bursts lie in these planes we estimated Ei^o and Liso us- 
ing the spectral parameters reported in the literature and listed in 
Table [1] excluding the two GRBs fitted with a single power law 
(i.e. with an unconstrained ^^^i^). Fig.[T]shows that the position 
of Fermi bursts is consistent with both correlations. We have fit- 
ted the 100 pie-Fermi GRBs, the 10 Fermi GRBs and the com- 
bined sample of 110 GRBs with the least square method. The 
best fit spectral parameters (normalization K and slope 6) and 
the probability that the correlation is by chance (P) are reported 
in Table |2] 



4. Time resolved spectra: data analysis 

The data of GRBs detected by the GBM since August 2008 are 
pubHcly available at http://fermi.gsfc.nasa.gov/ The 
GBM consists of 12 Nal and 2 BGO scintillation detectors differ- 
ently oriented so to derive the GRB position through the compar- 
ison of the count rates of the different detectors. The Nal cover 
the low energy spectral domain from 8 keV to ~1 MeV while 
the BGO detector (a factor 10 thicker than the Nal) are sensitive 
in the 0.2-30 MeV energy range. The GBM uses several trigger 
algorithms to detect a GRB. These are determined by different 
choices of the integration timescales and energy range. A trig- 
ger algorithm similar to that of BATSE Large Area Detectors 
(LAD) operates in the 50-300 keV band with a minimum signif- 
icance of 5 Act for the detection of the GRB over the background 
(which has typically a rate of 300-350 counts/s in this range - 
e.g. Meegan et al. 2009) on the 16 ms-4 s timescales. An event is 
flagged as a burst if at least two Nal detectors are simultaneously 
triggered. 

For our analysis we considered the TTE (Time Tagged 
Event) files containing the counts in 128 energy channels rel- 
ative to the burst period. We considered the TTE files of the two 
Nal detectors triggered by the GRB. From these we extracted 
the GRB light curve and the time resolved spectra with the gtbin 
tool (as part of the ScienceTools-v9r8p2-f ssc-20(59®225). 
Light curves were extracted by summing the count rates over 
the 8 keV-1 MeV energy range of the Nal detectors that were 
triggered and the 200 keV-30 MeV energy range of the 2 BGO 
detectors. Time bins of 1 s were adopted for all bursts. Light 
curves were rescaled in time to the trigger time of the GRB. 

Fig. |2] shows the burst light curves. In the case of GRB 
081007 the data present in the GBM catalogue 081007(224) re- 
fer to GRB 081007B which had very low statistics (Bissaldi et 
al. 2008). We could not find the data directory of GRB 081007 
which triggered Swift (Baumgartner et al. 2008) and Fermi 
(Bissaldi et al. 2008) 121 s after GRB 081007B. We only re- 
port in Tab.[T]the time integrated spectral results of GRB 081007 
(Bissaldi et al. 2008) . We note also that Fermi and Swift both 
detected a GRB on November 18 2008. These are, however, 
two events that do not coincide both spatially and temporally: 
RA=82.6° and dec=-43.3° is the location (with an uncertainty 
of 1.6' - Palmer et al. 2008) of the event detected by Swift at 
14:56:36 UT (Hoversten et al., 2008) while the Fermi event is 
located at RA=54° and dec=-50.4° (with an uncertainty of 2.9°) 
and was detected at 21:00:53.5 UT (Bhat et al., 2008a). Only 
for the Swift GRB the redshift is measured (z=2.58 D'Elia et al. 
2008, no E°\ measured). 

peak ^ 

For the scope of this work we do not subtract the background 
count rate to the light curves, they are only used here to show 
the time intervals within single GRBs selected for the extrac- 
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Fig. 1. The lOFermi-GRBs with redshift and E°ll^^ in the fipeak- 
Eiso (top panel) and Speak - hso (bottom panel) planes. Their 
position is compared with the correlations defined by the sample 
of 100 GRBs detected by other instruments (grey filled circles), 
i.e. pie-Fermi sample. The solid line is the best fit to the pre- 
Fermi sample, while the dashed lines represent its 3cr scatter 
The dot-dashed line is the best fit to the 10 ferm/-GRBs. Fermi 
bursts appear to be fully consistent both with the Speak - Eiso and 
the Speak - Siso correlations. 



tion of time-resolved spectra. For 3 GRBs (080905, 080928 and 
090423) the count rate is too small to perform a time resolved 
spectral analysis. In all the other cases (except GRB 081007) we 
could divide the light curve into time intervals as indicated by 
the dotted vertical lines of Fig. |2l For GRB 080916C we used 
the time resolved spectral analysis results reported by Abdo et 
al. (2009) who combined the GBM and the LAT data. Spectra 
of the triggered detectors were extracted over the selected time 
intervals defined with the gtbindef tool. The spectrum of the 
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Fig. 2. Fermi light curves of the bursts with known redshift detected by the GBM from August 2008 to July 2009. Light curves have 
time resolution of 1 s and are not background subtracted. The vertical dotted lines mark the time intervals selected for the extraction 
of the time resolved spectra. The light curve of GRB 081007 is not reported because the data of this burst could not be found. 



background was extracted in a time interval after the burst in 
order to limit the GRB contamination. Rebinning (with the grp- 
pha(v3.0.1 ) tool) was applied to each spectrum in order to have 
a minimum of 40 counts per energy channel. 

The response files corresponding to each detector were used 
for the spectral fitting. Spectra were analysed with Xspec(vl2) 
in the range 8 keV to ~1 MeV. For most GRBs the spectra 



of the two (or more) Nal triggered detectors were jointly fit- 
ted with a cutofl'-power law model (CPL) of the form F{E) = 
E" exp(-E/Ef)) with a free normalization constant for the spec- 
tra of the two or more detectors jointly fitted. Most of the fits 
give for this constant a value close to 1. The CPL has been 
widely used to fit the spectra of GRBs and in particular the 
time resolved spectra (Preece et al. 2000; Ghirlanda, Celotti & 
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Fig. 3. Time resolved fipeak - hso correlation of Fermi GRBs with known redshift. Small size symbols show the evolution of fipeak 
vs Liso- Large size symbols are the location in the fipeak - Lino plane of the corresponding bursts when the time integrated spectra 
(Tab. [TJ are considered. For GRB 081007 and GRB 090423 only the time integrated spectrum is available. The solid and dotted 
lines represent the iSpeak - ^^iso correlation and its 3cr scatter, respectively, as obtained with the pre-Fenni GRBs (Ghirlanda et al. 
2009). The dot-dashed line represents the fit to the 10 Fermi GRBs (time integrated) and the triple-dot-dashed line is the fit to the 
5 1 time resolved spectra. 



Ghisellini, 2002; Kaneko et al. 2006). We only analysed the 
spectra of the Nal detectors because in most bursts the inspec- 
tion of the BGO light curve did not show any evident signal. 
Morever, the lower sensitivity of the BGO detectors with respect 
to the Nal ones would lead, in extracting simultaneous spectra 
for a joint fit, to a much smaller number of time resolved spec- 
tra. 

For 8 GRBs of Tab. [T] we performed a time resolved spectral 
analysis. Note that also the BATSE time resolved spectra are of- 



ten fitted with the CPL function (e.g. Kaneko et al. 2006; Nava et 
al. 2008). For the purpose of comparing the spectral evolution of 
different bursts in the fipeak - i-iso plane the use of the same spec- 
tral model ensures that the possible biases, e.g. the overestimate 
of the peak energy with respect to the Band model (Band et al. 
1993), is a common systematic effect of all the analyzed spectra 
(e.g. see Kaneko et al. 2006). In any case, we also verified if our 
time resolved spectra could be consistent with the Band model 
finding that, in most cases, we could not constrain the high en- 
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Fig. 4. Left panel: the rest frame energy £0 - £'peak/(2 + a) as a function of luminosity for the time resolved spectra and for 
the corresponding time integrated spectra (same symbols as in Fig. O. Right panel: the low energy spectral index a versus the 
luminosity Liso- 



ergy power law slope of this model. This also motivated us to to 
choose the minimal simplest model, i.e. the CPL. 

We computed the isotropic luminosity of each time resolved 
spectrum by integrating the best fit spectral shape over the rest 
frame 1 keV-10 MeV energy range. Table [3] reports the results 
of the time resolved analysis: Col. 2 and Col. 3 give the start 
and stop times of the time resolved spectra. Col. 4 and Col. 5 
give the photon spectral index a and the characteristic energy Eq 
(with their 90% significance errors) respectively. 



5. Results 

5.1. Evolutionary tracks 

Fig. [3] shows the evolutionary tracks of the 8 Fermi GRBs with 
redshift for which time resolved spectral analysis was possible. 
The number of time resolved spectra extracted per burst depends 
on its total fluence. Our guideline in defining the time inter- 
vals was a trade-off between the need to follow the rise and 
decay phases of the single pulses within the light curve and to 
have enough signal in each time resolved spectrum to constrain 
its spectral parameters. In two bursts (GRB 090424 and GRB 
090618) a more dense sampling of the light curve is possible 
given their large fluence. Here we have limited the number of 
time bins to have a more neat evolutionary track in the fipeak-i'iso 
plane, while in the next sub-section we will discuss their time 
evolution in full detail. For GRB 090323 the time resolved spec- 
tra up to ~ 100 s after the trigger are best fitted by a simple 
power law model, while only after 100 s a curved model (CPL) 
can be fitted and the value of the peak energy can be constrained. 
Therefore, for GRB 090323 we show in Fig.[3]only the spectral 
evolution of the final peaks (those after 100 s in the correspond- 
ing light curve of Fig. |2]l. 

Fig. |3]shows that the prompt spectrum evolves in a well de- 
fined way, and that fipeak and Liso are correlated. The entire evo- 
lutionary tracks of all the 8 bursts we studied lie within the 3cr 
stripe of the scatter of the iSpeak - i^iso correlation defined by 
the time integrated spectra of different bursts. Fig.|3]shows data 
points that are associated both to the rising or descending part 
of pulses, although they are difficult to distinguish. In any case, 
there appears to be no difference between the rising and decay- 



ing part of the pulses, and this important point will be discussed 
in more detail below. 

Since fipeak is a derived quantity, being - fipeak = £^o(2 + a), 
it is interesting to verify if the E' , - L' correlation is the re- 

^ J peak ISO 

suit of underlying correlations of Eq and/or a with Liso. The left 
panel of Fig. |4] shows the rest frame energy £0 - £^peak/(2 + a) 
as a function of the luminosity for the time resolved spectra and 
for the corresponding time integrated spectra (same symbols as 
in Fig. [3]). The right panel of Fig.|4]shows the low energy spec- 
tral index a versus the luminosity Liso. We can see that a shows 
no correlation with Liso when considering the ensemble of bursts 
(but see below the case of the single pulse of GRB 090618), 
while Lo does. We can conclude that it is indeed Lpeak (or Lo, 
i.e. the e-folding energy of the CPL model) the spectral param- 
eter correlating with Liso. From the right panel of Fig. |4] we also 
see that the low energy photon index of the time integrated and 
time resolved spectra of the analyzed Fermi GRBs are consis- 
tent with the -2/3 limit predicted by synchrotron emission while 
they violate the -3/2 synchrotron limit for fast cooling electrons 
(Ghisellini, Celotti & Lazzati 2000). 

We analyzed the E* , - L' correlations obtained with the 

^ peak ISO 

time resolved spectra. In Table|2]we report the normalization and 
slope (K and 6, respectively) and the chance probability P of the 
correlation. In total we have 5 1 time resolved spectra for the 8 
GRBs. They define a correlation in the Lpeak - ^iso plane with 
a slope 0.37+0.05 (triple-dot-dashed line in Fig. [3]) consistent 
with the slope of the correlation defined with the time integrated 
spectra (i.e. 0.39+0.03 for the 100 pre-Ferm/ GRBs - solid line 
in FiglU- or 0.40+0.02 for the total of 1 10 GRBs including the 
Fermi events - dot-dashed line in Fig.[3]l. 



5.2. Rising and decaying ptiases of single pulses 

GRB 090424 consists of two peaks overlapped with a total dura- 
tion of 6 seconds. GRB 090618 has a complex light curve made 
of a smooth precursor followed by two intense peaks partially 
overlapped (see Fig. |2| for a total duration of 150 s. Their high 
count rates allow a time resolved spectral analysis with a dense 
time sampling. These two GRBs are well suited to study how 
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090424 



090618 




Luminosity [erg/s] 




Luminosity [erg/s] 



Fig. 5. High time resolution spectral analysis of GRB 090424 (left) and of GRB 090618 (right). Top: light curve (0.1s resolution) 
and time intervals of the extracted spectra (vertical dotted lines). The coloured symbols represent the spectra highlighted in the 
panels below. Spectral evolution panels (from top to bottom): rest frame peak energy Zipeak, rest frame characteristic energy Eo and 
photon index a versus isotropic luminosity. The filled circles and the open squares identify the first and the second peak of the light 
curves. The colour code marks the rise/decay phase within each peak. 



the spectrum evolves during the rise and decay phases of their 
pulses. 

We divided the time interval of the duration of GRB 090424 
so that each extracted spectrum had a signal-to-noise ratio 
S/N>30, integrated over the 8 keV-1 MeV energy range. This 
gives a total of 42 spectra distributed in the ~6 s of duration 
of GRB 090424. For GRB 090618 we required a S/N ratio of 
50 and obtained 34 time resolved spectra in the 150 s of its du- 
ration. These spectra were extracted and analysed as described 
in §3. In Fig. |5]the light curve (with 0.1s time resolution) and 
the three panels of the correlation between the spectral param- 
eters (iipeak, £^0 and a) and the luminosity Liso are shown for 
GRB 090424 (left) and GRB 090618 (right). In these plots we 
have marked with different symbols the spectral evolution of the 
different peaks identified in the light curve and with different 
colours the rise and decay phase of single pulses. The different 
peaks have also been fitted individually in the fipeak - i'iso plane 
(second panels of Fig. |5]l and the results are shown in Table |2] 



We note that these 4 peaks (two of GRB 090424 and two of 
GRB 090618) define a very tight E^^k - Li^o correlation with 
slope between 0.45 and 0.6. 

These plots indicate that the rising and decaying phases are 
indistinguishable in a iipeak - ^^iso plot. Intriguingly, we also find 
that for GRB 090618 there is a correlation between a and Liso. 
Since also Eq correlates with Liso (as for the other bursts), this 
results in a even tighter E'^^^^ - L[^^ correlation for the two pulses 
of this GRB. 

6. Discussion & Conclusions 

The largest sample of long GRBs with measured redshift and 
^peak collected recently (e.g. Ghirlanda et al. 2009; Nava et 

al. 2008) defines a strong correlation Lpeak °^ L^^^^ with scatter 
cr ^ 0.26. A similar strong correlation exists between the peak 
energy and the isotropic energy, i.e. Lpeak^^ E^^j^l- The time inte- 
grated spectra of the 8 Fermi GRBs with measured redshift (open 
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symbols in Fig. [T]) are consistent with both the fipeak - £^iso and 
the iipeak - iiso Correlation defined by 1 00 pre-Fe rm/ bursts (grey 
filled circles in Fig.[T]i. In the fipeak -i'iso correlation, Zipeak is that 
of the time integrated spectrum and Liso is computed (as done by 
Yonetoku et al. 2004) using the value of iSpeak of the time inte- 
grated spectrum and the peak flux of the GRB. This luminosity is 
somehow different from the luminosity of the brightest phase of 
the burst, i.e. the luminosity of the peak spectrum. A discussion 
about these two definitions of Liso can be found in Ghirlanda et 
al. (2005). 



6. 1 . Selection effects 

The sample of bursts defining the ZSpeak - £^i,so and iSpeak - ^iso 
correlations is heterogeneous: Nava et al. (2008) considered 
83 pie-Fermi GRBs detected by different instruments (BATSE, 
BeppoSAX, Hete-II, Konus-Wind and Swift) since 1997, a more 
recent update (Ghirlanda et al. 2009) consider 100 GRBs. Here 
we have added the 8 Ferm/-GRBs with measured redshift. 
Adding new GRBs to the above correlations does not represent 
a secure test of their physical nature especially when there is the 
suspect that these correlations are the result of instrumental se- 
lection effects. Different instrumental selection effects could be 
biasing the samples of bursts used to define the ZSpeak - ^iso and 
£peak-iiso Correlations. Butler et al. (2007) (see also Shahmoradi 
& Nemiroff 2009) argued that the spectral-energy correlations, 
in particular the Lpeak - ^iso and the Lpeak - Li&o defined by the 
time integrated spectra, are the effect of the trigger threshold, 
therefore having no physical relevance for the understanding of 
the GRB emission mechanism. In Ghirlanda et al. (2008) we in- 
vestigated this issue by studying instrumental selection effects 
possibly biasing a sample of 76 bursts (updated to Sep. 2007) 
with measured redshifts. This sample defines strong correlations 
in the observer frame E°^^, -F and E°'°^, -P planes (F and P are 

peak peak ^ ^ 

the fluence and peak flux) where the instrumental selection ef- 
fects can be studied. Two selection effects were considered for 
the detectors on-board BATSE, BeppoSAiX and Swift: the trig- 
ger threshold, i.e. the minimum flux a burst must have to trigger 
a given detector, and the spectral threshold, i.e. the minimum flu- 
ence a burst must have in order to constrain its spectral param- 
eters (in particular the peak energy £pgj||j). Our results indicate 

that: (i) both the selection effects are functions of E^^^ but the 
spectral threshold is dominating over the trigger threshold, im- 
plying that the Lpeak - £^i,so or Lpeak - ^-iso Correlations are not due 
to the trigger threshold; (ii) the Swift spectral threshold is biasing 
the Swift GRB sample with redshift added to the Lpeak - ^iso cor- 
relation in the last three years; (iii) selection effects are present 
but they do not determine the spectral-energy correlations. 

Another way to test the incidence of selection effects on the 
Epeak - Eiso and Lpeak " i'iso Correlation is to verify how much the 
slope, normalization and scatter change by separating the het- 
erogenous sample of GRBs with redshift into sub-samples of 
bursts detected by different instrument (e.g. Butler et al. 2007). 
Butler et al. 2007 found similar slopes but different normaliza- 
tions of the Lpeak - £^iso Correlation by considering the pie-Swift 
and the Swift sample (but see Amati et al. 2009). Unfortunately, 
Swift bursts have a very narrow range of limiting the ro- 
bustness of this test (Nava et al. 2008). Furthermore, different 
instruments like BATSE and BeppoSAX can have very similar 
detector thresholds and bias, in a similar way, the GRB samples 
that they detect. 



6.2. Time resolved E' , - L' correlation 

peak ISO 

All the spectral-energy correlations have been derived consider- 
ing the time integrated GRB spectral properties. By studying the 
spectral evolution of the 8 Fermi GRBs with measured redshift 
we also find that a correlation £' ,-L' , between the rest frame 

peak ISO' 

peak energy Lpeak and the bolometric isotropic luminosity, exists 
within individual bursts (Fig. [3]). This E'^^^^^ - U^^ correlation can 
also extend over two orders of magnitude in both iSpeak and Liso 
within the duration of a burst. The evolutionary tracks defined 
by the 8 Fermi GRBs lie in the upper part of the Lpeak - i-iso 
correlation. This could be due to a systematic underestimate of 
the luminosity in time resolved spectra with respect to time in- 
tegrated spectra (which are used to define the Lpeak - ^iso corre- 
lation plotted as a solid line in Fig.|3]l. Indeed, the time resolved 
spectra are more frequently fitted with a CPL model which lacks 
the high energy power law component of time integrated spectra. 
The finding of a £' ,-L' correlation within individual 

^ peak ISO 

GRBs, consistent with the Lpeak -i'iso correlation defined by time 
integrated spectra, is the strongest argument in favour of a phys- 
ical origin of this correlation and the strongest argument against 
instrumental selection effects biasing the observed correlations. 

6.3. Interpretations of the spectral-energy correlations 

A convincing way to ensure the reality of the spectral-energy 
correlations would be to find a robust physical interpretation. 
The proposed interpretations of the Lpeak - ^iso and Lpeak - £^iso 
correlations can be divided into two classes: (a) kinematic in- 
terpretations in which the link between Lpeak and Liso is estab- 
lished by the configuration of the emission region, i.e. a uni- 
form jet observed at different angles (Yamazaki et al. 2004), 
an inhomogeneous jet model (e.g. Nakamura 2000; Kumar & 
Piran 2000) made up of multiple sub-jets or emission patches 
(Toma et al. 2005) or a ring-shaped emission region (Eichler & 
Levinson 2004); (b) radiative interpretations in which it is the 
emission mechanism of the prompt phase to link Lpeak and Liso 
as in the case of a spectrum dominated by a thermal component 
(Meszaros & Rees 2007; Ryde et al. 2006; Thompson, Meszaros 
& Rees 2007), in the case of photospheric emission dominated 
by magnetic reconnection (Giannios & Spruit, 2007) or when 
the emission is synchrotron radiation from the external shock 
(Panaitescu et al. 2009). 

The common feature of the kinematic models in reproducing 
the Lpeak - i'iso or Lpeak " ^iso Correlation is the viewing angle 
under which different GRBs are observed. Both the off-axis and 
the sub-jet models need to assume the existence of a on-axis 
correlation between the peak energy and the luminosity, whose 
origin could be instead related to the radiative process. Indeed, 
the kinematical models that (under some assumptions about the 
typical jet opening angle distribution) succeed in reproducing 
£^peak<^ i'fso ^^ould Still explain a similar correlation within indi- 
vidual GRBs, i.e. a time dependent correlation E' , -L' , which 

' ^ pea k ISO 

can extend over 2 orders of magnitude (e.g. Fig.|4|. 

The simplest way to explain the Lpeak - i^iso correlation is 
to assume that only the bulk Lorentz factor F changes. Since 
^^peak'^ r and L oc F^, we recover Lpeak*^ L^^^. But this assumes 
that, in the comoving frame, both Lpeak' and L' are the same even 
if different F-factors are required, and this seems unlikely (both 
when considering the Lpeak - i-iso correlation defined by different 
bursts or the E' ^ - L' correlation holding within individual 

peak ISO ^ 

GRBs). 
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If the emission is due to the synchrotron process, the peak 
frequency iSpeak'^ ^r-^p (where is the random Lorentz factor 
of the electrons emitting at the peak) and L oc NB^T^y^, where 
is the number of the electrons having y^. Therefore a change 
of the quantity TB, maintaining the same y^ and A^, would give 
iipeak'^ L}^^- But the prompt emission almost surely occurs in 
the fast cooling regime, implying that the resulting synchrotron 
spectrum cannot be harder than L{E) oc £ '^^ (Ghisellini et 
al. 2000) while we observe (e.g. Preece et al. 2000; Ghirlanda, 
Celotti & Ghisellini 2003) also in the Fermi GRBs (Fig.Hright 
panel) much harder spectra. Furthermore, it is seems hard to 
maintain the same and y-p while changing TB. 

Quasi-thermal Comptonization could well explain the fact, 
found in GRB 090618, that both and a coiTelate with L. 
In fact, if the seed photons for Compton scattering remain the 
same, an increase of the plasma temperature would increase the 
Comptonization parameter y ~ 4r/rre/(meC^), producing both 
an harder spectrum and a larger fipeak (j is the optical depth 
of scattering electrons, and Tg their temperature). On the other 
hand, for likely bulk Lorentz factors V ~10~-10'', the comoving 
temperature is below 1 keV, implying r > 10^ to reach the re- 
quired y ~ 10, needed to account for the observed flat spectra. 
With this values of t the resulting spectrum would saturate to 
a Wien-like spectrum, not to a cutoff power law. A very large 
value of T would also lengthen any variability timescale. 

It has been suggested (e.g. Borgonovo & Ryde 2001; Ryde 
& Petrosian 2002) that the off-latitude emission that follows an 
abrupt switch-off of the fireball introduces a spectral-energy 
dependence, since the observer sees progressively less beamed 
(and less blue-shifted) emission. However, this could explain 
only the decaying phase of the pulse. In this paper we analyzed 
the spectral evolution of two of the most intense bursts in our 
sample: they allow to make a dense sampling of their light curves 
in order to extract time resolved spectra. These are GRB 090424 
and GRB 090618 (Fig.|4]right and left panel, respectively). Our 
findings indicate clearly that there is the same E' . - L' cor- 

^ J peak ISO 

relation during the rise and the decay phase of different pulses 
within these two GRBs. The time evolution is so that during the 
rise phase both fipeak and Liso increase to the maximum value 
and during the decay they decrease along the same evolution- 
ary track that they followed during the rise phase. In the case of 
GRB 090618 these considerations are valid also for the correla- 
tion between a and Liso (bottom right panel of Fig.|4]i. 

In the attempt to explain the fipeak - £^iso correlation, 
Thompson Meszaros & Rees (2007) pointed out the importance 
of shear layers shocks to extract a large fraction of the bulk ki- 
netic energy of the fireball, leading to a black body spectrum. 
The same arguments could be used to explain the fipeak - i^iso 
correlation, in diff'erent bursts and the £' , - L'. correlation 

' peak ISO 

within individual bursts as well. The problem with this inter- 
pretation is that one of the key assumption of their scenario is 
that the value of the bulk Lorentz factor in the dissipation region 
must be fine tuned (it must be of the order of l/0j, where flj is 
the opening angle of the fireball). This assumption is relaxed in 
the "reborn fireball" scenario (Ghisellini et al. 2007), but there 
remains to be explained why so few bursts have pure black body 
spectra (Ghirlanda, et al. 2003), and, even when adding a power 
law component (Ryde et al. 2005), its slope is too soft to explain 
low energy data (in the keV band), as shown by Ghirlanda et al. 
(2007). 

We can conclude that new ideas are called for explaining 
what emerges to be a general and well defined property of the 
prompt emission of GRBs. 
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-1 75+"-2 

1 -5+0.11 
"'^•-'-0.12 

-1 03+"' 
-81 

-0.96t° 

" -p.os 
-0 97+" 12 

-1 13+"' 

-1 04+"' 

-1-38S 



106+Jj 
88.22^7 

59+23 

30+^^ 

-'"-10 

302+^^1' 
176!«^ 
130!j^ 
92+^5 

1^7+498 



50+^f 

119+81 

145!^^ 
158+1 

124+^^ 

161+^^^ 

748+g2 
838+lf 

581!j« 
7741000 

369 
441+227 
^^^-126 

448+?85 



103!^« 



122(131) 


0.42 


156(148) 


0.4 


164(175) 


0.17 


225(209) 


0.06 


73(70) 


0.54 


55(59) 


1.39 


68(60) 


1.29 


61(50) 


0.7 


82(89) 


0.15 


105(88) 


0.64 


103(91) 


2.84 


108(90) 


5.0 


108(88) 


5.36 


96(91) 


1.83 


123(89) 


1.0 


87(88) 


2.11 


84(90) 


5.8 


118(87) 


4.9 


127(88) 


6.2 


103(90) 


9.8 


129(90) 


14.7 


93(87) 


3.4 


105(91) 


0.83 


144(91) 


0.66 


98(90) 


0.61 


98(87) 


0.27 


102(87) 


0.36 


129(86) 


0.38 


114(86) 


0.43 


99(87) 


0.38 


106(87) 


0.05 


135(84) 


0.9 


168(181) 


1.6 


215(166) 


1.45 


167(167) 


3.0 


91(102) 


1.94 


180(154) 


1.52 


118(119) 


2.35 


144(132) 


1.54 


169(159) 


0.25 
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Table 4. Spectral evolution. Continued 



GRB 


tl 


t2 


or 


Eo 


/3 x'idof) 


F-6 




s 


s 




keV 




erg/cm^ s 



090424(592) 



090618(353) 



-1 
1 


1 
2 


-0.81!«;| 
-0.83!°| 


2 
3 
4 


3 
4 
6 


rv Q1+O.0S 

-Q09 
-0 79+"' 11 

-0 n+o-o'' 


6 


10 


-^•ii_0.4 




3 

14 

50 

60 

63 

67 

70 

75 

80 

85 

88 

100 

114 



3 

14 

40 

60 

63 

67 

70 

75 

80 

85 

88 

100 

114 

130 



.1 ig+o.n 

-i.U8_oo7 
1 4^+0.08 
-0.08 
1 Q+0.06 
"^•-'-0.06 
-1 06+'''"' 

.1 0+0* 

^- -0.03 

07+«.04 

"l--l--'-0.04 

-1 21+"'-'''' 
1 io+8:8^ 

"^•^-'-0.06 

-1.26+! 
-1.5+"- 
-1.6 



:8f 



103!|« 

159+1 

105!^'' 
137+'^ 

•^^-18 

612+^"! 
258!f 
153!^^ 
247+53 

^^'-40 

289^5^ 
319 
237_ 
172^1^ 
125+'^ 
170+j^ 

130+!f 
8i:l^ 



80(70) 


3.4 


71(73) 


11.4 


70(66) 


6.5 


84(59) 


4.4 


104(76) 


5.78 


62(62) 


0.4 


124(120) 


2.1 


296(188) 


1.64 


267(214) 


0.76 


246(187) 


1.88 


166(150) 


5.12 


296(183) 


11.7 


204(167) 


8.9 


219(167) 


5.0 


280(256) 


2.76 


203(167) 


5.0 


153(144) 


5.0 


219(188) 


2.1 


233(189) 


1.35 


247(193) 


1.68 



